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>v ABSTRACT 

s. 

q . The plasma theory of electrostatic waves generated collisionlessly when 

neutrinos pass through degenerate electrons is developed. It is consistent with 

quantum statistical field theory in the degenerate limit. The reactive case 

can also be treated. For neutrino distributions considered herein, neither case 
OO , ' 

\Q . develops significantly during e.g. the neutrino burst in a supernova explosion. 

(N 

in 
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Q\ . 1. Introduction 

Oh A number of recent papers have speculated on the collisionless heating of thermal 

plasma behind a Type II supernova shock by the intense flux of neutrinos emitted from 
the supernova core. The idea is attractive since the total energy radiated by a supernova 
in neutrinos (and antineutrinos, ~ 10 53 — 10 54 erg) far exceeds that in the thermal plasma 
(~ 10 51 ergs), and so a small transfer of energy from neutrinos to electrons and nuclei, in 
excess of that already accounted for by collisional processes alone, could go a long way to 
solving what has become known as the Type II supernova problem, i.e., that in models, the 
supernova shock following the core collapse and bounce loses energy to nuclear dissociation, 
stalls and eventually collapses back. Such an event would produce neutrinos, but not the 
optical emission for which supernovae are famed. 

Bingham and coworkers [[TJ derived a growth rate for the neutrino analog of the 
reactive two stream instability in hot (~ 3 x 10 9 K) plasma, treating the coupling between 
the neutrinos and the plasma by an effective potential ^/2Gpn e , following [§, where Gp 
is the Fermi weak interaction constant and n e is the electron density. A number of the 
assumptions in this paper were criticized by Hardy & Melrose, who used an approach based 
on quantum statistical field theory to derive a growth rate for the kinetic instability. Aside 
from rather elementary assumptions like the use by [[|] of a monoenergetic neutrino beam 
rather than a Fermi-Dirac distribution, and the treatment of the thermal plasma with a 
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cold electron dielectric tensor Q show that their result for the growth rate is lower than 
that of |T[ taken in the kinetic regime by a factor (l — u 2 / |k| 2 c 2 ) , where u and k are 
the wave frequency and wavenumber respectively, and c is the speed of light. This is an 
important point, since for plasmas where the electrons are relativistic or quasi-relativistic, 
to ~ |k| c, and this extra factor gives a big reduction in the calculated growth rate. Below 
we will present a plasma theory treatment of the instability that reproduces the quantum 
statistical field theory result in the kinetic limit, and gives a result similar to in the 
reactive limit. 

Another important caveat to these works (also emphasised by pi) is the collisional 
damping as electrons within the wave collide with ions. Even the most optimistic calculation 
gives a neutrino growth rate that is orders of magnitude lower than the collisional damping 
rates. Assuming that the wave mode actually exists, neutrino induced oscillations will be 
rapidly damped giving rise to plasma heating. In [p], ||] the plasma was assumed to be 
hydrogenic, at an electron temperature T e = 3 x 10 9 K and density n e = 10 30 cm -3 , giving of 
order A ~ 100 electrons per Debye sphere. However heavy elements are often the dominant 
component of supernova plasma. For example Cas A is thought to be mainly O 0,|. In 
a heavy element plasma A scales as Z~ 3 ^ 2 where Z is the ion charge. For fully ionized O, 
A ~ 10, and the so-called plasma approximation may be expected to break down. In this 
case no collective motions will exist, and no collisionless energy deposition can occur. 

For this reason this paper will concentrate on degenerate electron plasma found at the 
supernova core. At densities 10 34 — 10 38 cm -3 the electrons are highly degenerate (Fermi 
energies 10 - 300 MeV compared with kT e ~ 0.3 MeV) and collisional damping is strongly 
suppressed. The anisotropy in the neutrino distribution will be taken to to come from 
the natal kick that a pulsar may acquire during the supernova event. Various authors 
( || [7| and references therein) have considered the role of neutrino scattering from partially 
polarized neutrons in a 10 15 G magnetic field that may be present as an explanation for 
this kick, (due to parity violation, the neutrinos preferentially scatter to directions along 
the magnetic field). However a recent paper (which also gives a convenient review of the 
observations of pulsar kicks) || suggests that this is unlikely to be the actual mechanism. 
However it is clear that if a pulsar is born moving with a particular velocity with respect 
to the supernova remnant, no matter what the precise mechanism for this is, then the 
neutrino emission must be anisotropic. This comes about simply because the momentum 
carried off by the visible supernova/supernova remnant is generally much smaller than that 
of the pulsar, and consequently cannot conserve momentum by itself. Such a scenario might 
plausibly result from an off-centre detonation. For a pulsar moving with velocity 1000 km 
sec -1 an anisotropy of order 1% |J must be present in the neutrino emission. 
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2. The Neutrino-Plasma Interaction 

The work of treated the neutrino plasma interaction by an effective potential 
v^2GVn e , following Bethe [[J, who was considering the neutrino refractive index in the solar 
core and envelope, and hence was dealing with a non-relativistic plasma. For neutrino 
interactions with a relativistic plasma the effective potential should be generalized. The 
amplitudes Fw, Fz for neutrino-electron forward scattering for W and Z exchange are (in 



units with % = c = 1) M, 11 



1 w 



-G F _ 

—^ u vlt* I 1 - 75) u v u e ^ (1 - 75) u e 
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Uul^ (1 - 75) UyU^f (- - 2 sin 2 6 W - u e , (2) 



where u a , u a are final and initial wavefunctions for particles of species a, (m a their mass, 
Pai Pa,v their 4-momenta,) and sin 2 6w ~ 0.23 is the Weinberg mixing angle. A Fierz 
transformation has been performed in the expression for Fw- Only electron neutrinos 
interact via W exchange, whereas all three neutrino flavours interact with electrons through 
Z exchange. The axial vector terms do not contribute to the excitation of longitudinal 
waves, so Fw » Fz- Henceforward we consider only Fw-, putting 

Uulu (1 ~ 75) u u = 2u 1/ 7 (U « i/ = 2u u u u p u ^/m u (3) 

where the last term follows from the Gordon identity. For electrons the same treament 
gives u e u e (p^/m e — cr), where cr is a 3-vector composed of the Pauli spin matrices. 
Including wavefunction normalization factors, \Jm a / E a and and integrating over a electron 
distribution f e (p e ) (p a is the 3-momentum for particle species a) with no net polarization 
we get the following result for the effective potential for electron neutrinos 

V eff = V2G F J f e ( Pe ) (l - rf 3 Pe . (4) 

Obviously for non-relativistic plasma the Bethe result is recovered. Also for isotropic 
plasma, electrons with oppositely directed momenta will average to zero in equation (f|) 
recovering \p2Gpn e . However if a plasmon exists in the degenerate plasma, it must be 
carried by electrons moving at or very close to the Fermi velocity, and this scalar product 
will not necessarily average out to zero. 

Following the development and notation in [12| (pages 19-20) we estimate instability 
growth rates. The force F acting on an electron becomes 

F = -q (E + v e x B) - V2G F V J f v (p„) (l - ^jf) d 3 p, (5) 
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where E and B are electric and magnetic fields, uo and k are the frequency and wavevector 
of oscillating quantities, — q is the electron charge and v e the electron velocity. This gives 
an extra current density Jj due to the neutrino interaction 

J, = V~2 q G F jvju (p.) (i - ^) d ^^M^ VldW (6) 

Putting V/„ = VSf v = ik6f v , pjE a = v a and 

5f u (p„) = -V2G F ( 5f e ( Pe ) (1 - • v e ) d 3 p e k ' df : /dPu , (7) 

J LO — K ■ ~V U 

with J, = o"jj-Ej and ifjj = 1 + ioijjuj we arrive at an expression for the extra term in the 
dielectric tensor Kij due to the neutrino interaction; 

5K,= ^ J J J { ( "- k Ze^ } Ik (1 - ^ • V e ) dV„ 

X^(l - V, • V e ) rfV^f ^Pe- (8) 

Specializing to longitudinal waves, 

xt S^( 1 -^-v.)<< 3 P^i|i^fd 5 p s . (9) 

This expression is evaluated by noting that integrals of the form / v e k ^^{^ Pe ^ 3 Pe are 
parallel to k for a spherically symmetrial electron distribution function f e (p e ), and those 
/ v e k ' 8 y Pe (i 3 Pe = 0. Hence we may put 

v„ • v e -> Vj, • k v e • k/ |k| 2 (10) 

and perform the integration using the standard results in Table 1, putting 
<9/e/<9|p e | = —2/ (27r) 3 x <5 ( |p e | — Pf) and cj^ = n e q 2 'j^pm e . We have neglected the 
variation of the relativistic 7 inside the integral and have approximated it by its value at 
the Fermi surface. The full result is 

2|k| v\\ \k\v f \uj - |k| v F J J 

with 

A= 2G , ^/fl-^^) 2 |^f2-^logf^^ N ) N ) 

V |k| |k|y 2p F v F \ \k\v F to \u)-\k\v F ) ) 
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and the subscript F denoting values of variables at the Fermi surface. If we let A — > 
we recover the usual plasmon dispersion relation for degenerate electrons (cf This 
dispersion relation is plotted in Figure 1 for an electron density of 10 34 cm" 3 . The vertical 
dotted line indicates the wavevector where uj = uj cr it = |k|. To the right of this line uj < |k| 
and neutrinos may excite electrostatic waves by a kinetic instability. To the left of this 
line they are kinematically forbidden from doing so, and only the reactive instability may 
operate. This line moves to higher wavenumber (in units of the plasma frequency) with 
increasing electron density. 



2.1. The Kinetic Growth Rate 



To derive the instability growth rate in the kinetic regime we let uj — > uj + i5 with 

uj » 5 and put a) _^. v ► P ( - _^. v J — i^S (uj — k • v„), which gives, after taking imaginary 

parts, 



2nG 2 F u\k\ 2 



X 



-|k|V 



/k-f^-k-v^Kp,. 



(13) 
(14) 



This result is identical to that derived in H] by the methods of quantum statistical field 



theory, with the substitution of the appropriate residue factor for degenerate electrons [13 



I 1 da 

2uj duj 



- £J=£j([k|) 



2(uj 2 - 


k| 


2 v%) 




3uJ 2 - (uj 2 




k 





(15) 



where uj 2 = —a (uj) is the plasmon dispersion relation, and we have substituted from the 
dispersion relation to eliminate the logartithm. Similar results for spontaneous Cherenkov 
emission by neutrinos have been obtained by [14, 15|. It is apparent the the factor 
^1 — uj 2 / |k| 2 ^) by which this expression differs from the kinetic limit of the growth rate 
derived by fl|] is not related to electro- weak parity violation effects as implied by ||, 
but to a proper treatment of the relativistic plasma. In fact it is straightfoward to redo 



Table 1: Standard Integrals 

f 1 J^dx = -2 + f log (*±g) 

J ± Lo—kvx kv ° \ui— kv J 

rl j2^L dx = j£_ ( — 2 -)- 77- log 

J J- LQ—kvx kv \ kv y lv — kv J J 

•^1 lo—kvx^^ 3~^(fct>) ( ^ kv (cj— feu) J 

f 1 kvxi d x = <a. ( —1 4- ( ±l\ 2 (-2 + — log ( u)+kv ) )) 

Jl uj—kvx kv \ 3 \ kv ) \ kv & \w— kv J J J 
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the calculation of the neutrino-plasmon vertex in ]3| admitting right-handed as well as 
left-handed neutrinos (i.e. replacing u v j^ (1 — 75) u v with 2u v ^^u u ) The result is the same 
as for the pure left-handed neutrinos except for an extra factor of 2. 



2.2. The Reactive Growth Rate 

Turning to the reactive instability, considered in |]T[, where the pole in the integrand lies 
in a region where df u /dp u = 0. This will be the case in the region of phase space where the 
plasmon dispersion relation satisfies uj » |k|. Using the dispersion relation to eliminate 

the log inside the integral and putting 2 — jj^- log (^rj^fe) = — f — f ^4 F — • • • for 
uj » |k| elsewhere, we solve equations ([TT|) and ([12]) with K L = for to get 

- 2 — j(i + ! e S £ -+) 

x{l-^/(l-^p) 2 iSif<i 3 P„}=0 (16) 

This is similar to the expression derived in fll], but different by the factor ^1 — n^-p^r) in 
the integral. The integral in equation (|T6|) can be further simplified by taking the limit 
uj » |k|. The reactive growth rate is given by the complex roots (if any) to equation flllf) . 
Assuming that Jk • df u /dp u d 3 p u = (this is true for the distribution functions considered 
below), then to lowest order in |k| /uj a quadratic equation in uj with two real roots results, 
giving no reactive growth. Going to higher order a cubic equation in uj 2 results, again with 
no complex roots. A negative real root for uj 2 would give an imaginary uj, implying reactive 
growth at zero frequency. However such solutions should not be over-interpreted. At low 
frequencies electrons will couple to protons, which in turn will couple to neutrons by the 
strong interaction, and these couplings would need to be included in the dispersion relation 
to get a physical result. 



3. The Neutrino Distribution Function and Numerical Estimates 

To get numerical estimates of the magnitudes of the various growth rates we will 
consider a neutrino distribution function of the form 

f» (P") = n 17WTXI + ^nCOS0„) . (17) 

(2vr) exp (|p„| /k B T v ) + 1 

which comes from a thermal Fermi-Dirac distribution (the first term, with the neutrino 
chemical potential assumed zero) modified by scattering with an approximate differential 
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cross section 

d<7/cffi= <7 O (l + P n COS0) (18) 

with cr = GpEl/4ir 2 . The neutrino opacity, r = r E 2 ~ (i^/MeV) 2 cm 2 for typical 
neutron star parameters. The neutron polarization P n is given approximately by 

P n = 2x 1(T 5 (£/10 13 G) (3MeV/T n ) (19) 

where T n is the temperature of the neutron star matter. Consequently rP n « 1 for 
reasonable parameter values and equation ([11]) is justified. As remarked above this neutrino 
distribution function based on elastic scattering from partially polarized neutrons is unlikely 
to be the cause of pulsar kicks [[J . The reason for this is that the anisotropy in the neutrino 
distribution is wiped out by neutrino absorption by neutrons, unless a difference in the 
magnetic field of at least 10 16 G exists between the two opposite poles of the newly formed 
neutron star. However the neutrino distribution function that results when a pulsar receives 
a kick is likely to be of this form, so we proceed assuming equation fll?]) . For the plasma 
instabilities, the important feature is that neutrino scattering cross section rise with energy 
as E 2 . Thus the scattered (i.e. anisotropic) part of the neutrino distribution function will 
have a positive value of df v /dp v giving rise to postive kinetic growth rates in this region of 
phase space. 



3.1. The Kinetic Case 

The integral over the neutrino distribution function may be carried out putting 

k • df u /dp u = |k| cos pau/d \v u \ - |k| sin/5/ | P „| df v /de v (20) 

with cos/? = k • p u / |k| |p^|. For emission along 9 V = 0, (3 — —9 U and using the standard 
result / °° i™" 1 / (exp (x) + 1) dx = T (n) ( (n) (1 — 1/2™ -1 ) the integral evaluates to 



f 

dp/ v ~ " '" y ~ ^ 24 ' 240 



k • - — 6 {u - k • v„) d = — h — [k B T u ) T P n (21) 



for uj ~ |k|. The first term gives a negative contribution due to neutrino Landau damping 
by the thermal part of the distribution function. The second term gives rise to the 
instability. Thus 0.7tt 2 (k B T u ) 2 r P n > 1 is required. This evaluates to k B T u ~ 85 MeV, 
comparable to the gravitational binding energy per nucleon released in the collapse, for 
Pn ~ 2 x 10 -5 . At higher magnetic field strengths such as those thought to be present in 
magnetars |L6|, |17| (B ~ 10 15 G) the required neutrino temperature drops to 8.5 MeV. This 
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temperature (though probably not the magnetic field) is more consistent with neutrinos 
actually observed from SN 1987A JT8J. 

However as pointed out by |§ the kinetic growth rate will still be very small. 
The reason is that neutrinos may only excite waves with phase velocities in the range 
Vp — > 1. Reinstating the necessary factor of h~ 3 c~ 5 in equation (|T4]) and evaluating gives an 
approximate expression for the maximum growth rate 

5~1.8xl0~ 53 ^(l-^) 3 . (22) 

If 

for a neutrino temperature of 10 MeV, and the value of the integral in equation (|2~l"D taken 
to be + (/c^Tj,) 2 /24. This equation is evaluated for a variety of densities in Table 2. These 
values are far too small to be important during the few second neutrino burst during a 
supernova explosion. The lower densities in Table 2 may occur outside the neutrinosphere 
(the central region where neutrinos are trapped and are essentially isotropic, see j5]). In 
this case a beam instability of the type considered in [I], |3| would be more realistic. 



3.2. The Reactive Case 

With the same notation as above, the integral over the second term in equation (p0|) 
gives zero, while the first term evaluates to (hsT) 2 / 36. As expected this term gives no 
wave growth, since it represents the Landau damping of waves by the thermal neutrino 
distribution. We remark here that our analysis of the reactive instability is necessarily 
rather superficial. It would not be difficult to assume a neutrino distribution function that 
would give wave growth. However the question of whether such a distribution function 
is realistic for a supernova would still not be answered. For this reason the subsequent 
discussion will be concerned only with the kinetic instability. 



Table 2: Kinetic Instability Parameters 



n e (cm 3 ) 


1 — Vp/c 


7 


Up (rad/s) 


Merit (^p) 


5 (rad/s) 


10 32 


0.016 


5.6 


2.4 x 10 20 


2.11 


3 x 10- 7 


10 34 


7.5 x 10~ 4 


26 


1.1 x 10 21 


2.92 


3 x 10~ 9 


10 36 


3.2 x 10~ 5 


120 


5.2 x 10 21 


3.64 


3 x 10~ n 


10 38 


1.6 x 10~ 6 


550 


2.4 x 10 22 


4.21 


3 x 10~ 13 
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4. Discussion and Conclusions 

The fundamental reason that keeps the kinetic instability from developing significant 
wave growth is that the phase velocity of the waves supported by the plasma is very 
close to unity (the speed of light). Hence the factor of (l — oo 2 / |k| 2 ) in the growth rate 
becomes very small. The physical reason for this is that when the dispersive properties of 
the plasma approach those of the vacuum, longitudinal waves will cease to exist. With this 
in mind we will briefly consider the effect of the magnetic field on the plasmon dispersion 
relation. Again working from ref ]l2j (page 166, equation 10.21) we can derive an expression 
for K L with an ambient magnetic field. For simplicity we consider the case k • B = 0. 
For k||B, the problem is identical to the non-magnetic case discussed above. Putting 
dfe/dp e± = -2/ (2tt) 3 x 5 (p e± - ^P F -pt\\), J s -i (x) + J s+ i (x) = 2fJ s (x), and rewriting 
the sum over Bessel functions as — > Do° the expression is 

The electron cyclotron frequency is Q = q |B| /jpm e . For Q << |k| vp the sum over s can 
be replaced by an integral, whereupon its evaluation recovers the non-magnetic longitudinal 
dispersion relation. It is more interesting to consider the behaviour of individual harmonics 
over the range |k| vp << f2 — > |k| vp » fi. Then a behaviour something like the Bernstein 



modes discussed in [12] appears. For example for s = 1, as |k| — ► 0, uo — >• Qjjh where 



Vl\j H = Up + Q 2 is the upper hybrid frequency. However as |k| — > oo, u> — > Q ^ |k|, so a 
kinetic instability could have significantly higher growth rate than in the estimates above. 
Another possibility if Q > u p is a kinetic instability exciting electromagnetic waves, since 
for low frequencies Kp = |k| 2 jio 2 > 1. The discovery of pulsars with magnetic fields ~ 10 15 
G [ P| O opens up this possibility for degenerate electrons at densities 10 33 — 10 34 cm -3 . 
These novel scenarios will be investigated further in future papers. 

In conclusion then we have given a plasma theory for the so-called neutrino instability, 
that reproduces the results of quantum statistical field theory in the kinetic limit, and 
which gives more accurate expressions to use in the development of the reactive instability. 
The growth rates for longitudinal waves are insignificant in, e.g. the neutrino nurst during 
a supernova explosion, but we have been able to identify plausible conditions in which such 
wave growth might occur. 
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Fig. 1. — The dispersion relation for electrostatic waves in a degenerate electron plasma of 
density 10 34 cm -3 . The vertical dotted line indicates the wavevector where u = |k|, which 
is the onset of the kinetic instability regime. 



